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triangular lattice 



(N 
O 
O 
(N 

C 
3 



Ryotaro Arita, Kazuhiko Kuroki 1 , and Hideo Aoki 

Department of Physics, University of Tokyo, Hongo, Tokyo 113-0033, Japan 
1 Department of Applied Physics and Chemistry, University of Electro- Communications, Chofu, Tokyo 182-8585, Japan 

(February 1, 2008) 

We propose theoretically that a magnetic field can realize spin-triplet superconductivity in repul- 
sively interacting electron systems having strong ferromagnetic spin fluctuations. We confirm the 
general idea for the low-density Hubbard model on a triangular lattice, whose Fermi surface consists 
of disconnected pieces, by calculating the pairing susceptibility in a moderate magnetic field with 
the quantum Monte-Carlo method combined with the dynamical cluster approximation. 
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Triplet superconductivity is receiving current interests 
as one of the most fascinating new frontiers in super- 
conductivity. Although there are in fact mounting ex- 
perimental evidences for triplet pairing in various mate- 
rials including aii organic conductor TMTSHllB, heavy 
fermion systemsu, or a ruthenateuij, theoretical under- 
standing is far from straightforward. Among a vari- 
ety of possible approaches to realize triplet pairing, the 
mechanism which exploits the exchange of ferromag- 
netic spin fluctuation has long been studied, dating back 
to the pat&wiagnon interpretation of supepfluid 3 He in 
the 19ftl'sa0. Recent discoveries of UGe 2 B, ZrZn 2 B or 
UGeRhL3 where the superconductivity coexists with fer- 
romagnetism have invoked a renewed interest. 

One big question about the triplet pairing is how to 
raise its Tc- While singlet superconductivity mediated 
by antiferromagnetic spin fluctuation almost certainly 
enjoys its realization in high-Tc cuprates, triplet super- 
conductivity mediated by ferromagnetic spin fluctuation 
should in general have lower T c 's. This is mainly because 
spin a triplet pair can only exploit one component (e.g. 
longitudinal [zz] one for an S z — 1 pair) in the spin fluc- 
tuation as the pairing interaction, while a singlet pair 
can exploit all the three (two, transverse [-1 — ] as well as 
longitudinal) componentsOtl 

To overcome this difficulty, various-authors have pro- 
posed various ideas. Monthoux et al£3 have shown phe- 
nomenologically that the Ising-like anisotropy in (ferro- 
magnetic) spin fluctuationsj-Javors the triplet supercon- 
ductivity in the zz channcliil. Here we can notice that 
the anisotropy can be readily realized if we introduce 
an external magnetic field. Namely, when the ferromag- 
netic spin fluctuation is dominant, the spin susceptibility 
Xzz = J cIt(S z (t)S z (0)) is enhanced and an Ising-like 
anisotropy will be realized in a moderate magnetic field 
B (assumed here to be || z). Thus here we test the idea 
by studying the effect of a magnetic field to the triplet 
pairing instability in the ferromagnetic microscopic mod- 
els rather than a phenomenology. 

On the other hand, motivated h« the ferromamietic su- 
perconductivity in UGezIa, ZrZnJ and UGeRbH, Kirk- 



patrick et aZ.El have proposed that the ferromagnetic 
spin-fluctuation mediated superconductivity should have 
a higher T c in the ferromagnetic phase than in the param- 
agnetic phase. With field-theoretic techniques they have 
shown phenomenologically that the coupling of longitu- 
dinal spin fluctuations with transverse spin waves plays 
a crucial role to raise T c . 

Thus coexisting ferromagnetism and superconductivity 
should also be an interesting avenue, but theoretical real- 
ization of this is again not straightforward. Namely, the 
present authors have systematically studied the Hubbard 
model on various two-(2D) and three-dimensional(3D) 
lattice structures with the fluctuation exchange (FLEX) 
approximation, and have concluded that 2D systems 
are much., more favorable for superconductivity than 3D 
systemsEll. Physical reason for this has been identified 
in a much larger volume fraction of the phase space 
over which the interaction contributes to the Eliashberg 
equation. A similar conclusion has been obtained inde- 
pendently .by. Monthoux et al from a phenomenological 
calculatiorOlij. This should also apply to the ferromag- 
netic superconductivity, but a serious problem in 2D is 
the Curie temperature T c , which would upper-bound the 
ferromagnetic superconductivity, should be much lowerlill 
in 2D due to quantum fluctuations. 

This is exactly why we consider applying a magnetic 
field to an "almost ferromagnetic" 2D system to stabilize 
the ferromagnetism. So, in sharp contrast to a conven- 
tional wisdom that superconductivity and magnetic field 
are inimical to each other, we are here talking about a 
unique situation where the pairing arises due to the mag- 
netic fieldta. To repeat the idea, we consider a micro- 
scopic model close to ferromagnetism, which should be- 
come spin-polarized in a relatively small magnetic field 
even at high temperatures. Then we can possibly ex- 
ploit the favorable modes a la Kirkpatrick et al. without 
discarding the advantage of 2D. 

As for the microscopic model we take here the repulsive 
Hubbard model. While it has been shown, as mentioned 
above, that triplet pairing is ha*d to realize in the Hub- 
bard model in ordinary latticesO, a notable wayout has 
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been proposed by two of the present authorsta, where 
the triplet pairing is shown to become dominant if we 
consider systems having Fermi surfaces consisting of dis- 
connected pieces. There, we can insert the nodal linc(s) 
that are required to exist in a triplet gap function, in 
between the pieces, which greatly helps in enhancing Tq- 
A simplest example where such a Fermi surface occurs 
is the double-dipped band dispersion in the triangular 
lattice (FigJIJ) when the electron concentration is dilute. 
Ref.Ej has shown that an f-wave triplet pairing is indeed 
much favored than usual. However, it is still a subtle 
problem whether we can obtain a finite T c . Thus it is 
a theoretical challenge to clarify whether the triplet su- 
perconductivity in the zz channel (i.e., the spins of the 
pair || B with the d vector d _L z) can become dominant. 
If such a triplet instability is enhanced by applying an 
external magnetic field, it will open a new avenue for 
the triplet superconductivity^. We have calculated the 
triplet pairing susceptibility in magnetic fields, and show 
that the magnetic-field induced triplet superconductivity 
should exist, although we have not been able to quanti- 
tatively estimate Tq. 

From the motivation of the study, the method should 
be able to handle, in a microscopic way, higher-order 
processes such as the coupling between longitudinal spin 
fluctuations with transverse spin waves, and pairing fluc- 
tuations in low dimensions. Here we take the quan- 
tum Monte-Carlo (QMC) method in combination with 
the dynamical cluster approximation (DCA) formulated 
recently by Jarrell et allBlB. In DCA, the Hubbard 
model is mapped to a self-consistently embedded clus- 
ter of Anderson impurities rather than a single impu- 
rity considered, in the dynamical mean field approxima- 
tion (DMFA)E^I, so that DCA systematically incorporates 
nonlocal spatial fluctuations. Specifically, DCA+QMC, 
being a non-perturbative approach, automatically in- 
cludes the higher order processes. We have solved 
the the cluster problem generated by the DCA using 
the QMC methadcil with the algorithm proposed by 
Hirsch and Fyeo. We choose the cluster size N c = 
4x4, which is large enough to treat the anisotropic 
superconductivity^. As for the effect of the magnetic 
field, we concentrate here on the Zeeman effect by apply- 
ing the magnetic field B parallel to the plane to ignore 
the orbital effect. 

The Hubbard Hamiltonian is 

ija i ia 

in the standard notation, where tij is nonzero only be- 
tween nearest neighbors, h = g/J-sB is the Zeeman energy 
and the spin quantization axis taken to be || B. Here- 
after, the Coulomb repulsion U is set equal to the band 
width W, which is taken to be 2 here. 

The one-electron band dispersion, e — 2t[cos(k x ) + 
cos(fc x ) + cos(k x -I- k y )], for the triangular lattice (Fig. [I]) 
has two pockets for the Fermi surface when the band is 



less than quarter-filled (n < 0.5). For comparison with 
the case where ferromagnetic fluctuations dominate but 
the Fermi surface is simply connected, we also study the 
low-density Hubbard model on a square lattice with next- 
nearest transfers (t-t' lattice hereafter). Its dispersion is 
e = 2t[cos(k x ) + cos(fc x )] + 2t' cos(k x ) cos (k y ), for which 
we take t' = Q.5\t to make the the density of states di- 
verge at the, band-bottom to realize strong ferromagnetic 
fluctuationsllili23'L3. The Fermi surface for this lattice is 
a simply-connected star shape when dilute. To facilitate 
comparison with triangular lattice we adjust t to make 
W = 2(=U). 

Let us now present the results. We start by confirming 
that the systems considered here all have strong ferro- 
magnetic fluctuations. In Fig. ||, we plot the inverse 
spin susceptibility (x) versus temperature for the trian- 
gular and t-t' lattices. We can see that both systems ap- 
proach ferromagnetism (i.e., \ increases) for T — > 0. The 
tendency is stronger for a dilute filling (n = 0.2) with a 
disconnected Fermi surface than for a near-quarter filling 
(n = 0.4) for which the Fermi surface is almost simply- 
connected. 

In the presence of the Zeeman energy the triplet pair- 
ing susceptibilities, x^ ~ ( c j c { c T c t)> ~ ( c | c } c 1 c t)' 
X U ~ (cjcjcjcj.) with Sl otal = 1,0,-1 respectively, dif- 
fer from one another. Figure |^ plots these versus ex- 
ternal magnetic field at temperature T(= 1//3) = 1/32. 
For the two Fermi pockets in the triangular lattice we 
can consider an /-wave (with three nodes) gap function, 
<\> = x /3/32[cos(fc x /2) cos(fcp/2) sin((fe x - k y )/2)]. For 
this pairing symmetry all the nodes avoid intersecting 
the Fermi surface, so that the pair-scattering processes 
all around the Fermi surface(FS) has nonzero contribu- 
tions to the pairing vertex, 

Efc,fc' eFS ^(fc-fc , MfcMfc / ) 

which measures the tendency toward the pairing insta- 
bility. We can see that the /-wave susceptibilities of tri- 
angular lattice for n = 0.2 are indeed much larger than 
those for the t-t' lattice [with = ^/l/2 sin(fc x )]. x is 
larger for the dilute band (n — 0.2). These are consistent 
with the results of our previous calculation by FLEXEj. 

As for the magnetic field dependence, we can see 
that x is indeed enhanced in the presence of mod- 
erate magnetic fields in all the cases studied. The en- 
hancement, however saturates for a larger magnetic field 
(h ~ 1.0 x 10~ 2 ), so that the enhancement is peaked at 
some h(= /i ma x)- This is also the case with the most 
promising candidate, the low-density triangular lattice. 
This can be interpreted as follows. While the magnetic 
field does contribute to x from both the Kirkpatrick 
processes and the anisotropy in the spin fluctuations, the 
Fermi surface for the majority spin expands and becomes 
simply-connected when the Zeeman energy is too large. 
The nodes in 4> will then have to intersect the Fermi sur- 
face, which should degrade the pairings of majority spins. 



Therefore, the magnetic-field induced superconductivity 
can be expected only in the appropriate magnetic field. 

The "appropriate" magnetic field is dictated by the 
temperature. So we move on to the T-dependence of the 
pairing susceptibilities. We focus on x^{h = 1-0 x 10~ 2 ) 
and x(h = 0) in Fig. § x(h = 0) grows toward T = 
for the triangular lattice with n = CCi which is again 
consistent with the FLEX calculationta. As for the en- 
hancement of x (h = 1.0 x 1CP 2 ) above xQ 1 = 0), the 
enhancement decreases toward T = for the dilute tri- 
angular lattice. Extending the above argument we can 
see that this should be the case. For a fixed magnetic 
field the spin polarization becomes larger as the temper- 
ature is decreased, so that the majority-spin Fermi sur- 
face again becomes simply-connected. This implies that 
X^ continues to be enhanced at low temperatures only 
when we decrease the magnetic field accordingly. 

The situation is summarized in Fig. || which schemat- 
ically depicts the temperature-dependence of 1/x • A 
bunch of lines represent the behavior for various values 
of h, where x(Vx) a t each T first increases (decreases) 
and then reversed when h exceeds /i max , where /i max de- 
creases with T. If we assume that the derivative of x 
with respect to h continues to be positive around Tc at 
which x{h — 0) diverges, the bunch of lines forms an 
envelope in such a way that the Tc is enhanced for an 
intermediate value (h ~ 0.01) of the magnetic field. The 
value h = 0.01 corresponds to O(10T) for the band width 
- 0(1 eV). 

From the above argument a system in which the Fermi 
surface remains disconnected for larger Zeeman splitting 
may seem more favorable. However, we have to have 
an almost flat bottomed dispersion to realize ferromag- 
netism to begin with, so the persistent disconnectivity 
might be incompatible with this. Whether this can be 
overcome is an interesting future work. In conclusion, 
we have suggested an exotic possibility of a magnetic- 
field induced triplet superconductivity. This work was 
supported in part by a Grant-in-aid for scientific research 
and Special coordination funds from the ministry of edu- 
cation of Japan. Numerical calculations were performed 
on SR8000 in ISSP, University of Tokyo. 
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FIG. 1. Gray-scale plots for the band dispersion (left pan- 
els) and the dominant gap function (right) for the Hubbard 
model on the triangular lattice (a) and t-t' lattice (b). 
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FIG. 3. The triplet pairing susceptibilities versus magnetic 
field for the triangular and t-t' lattices for n — 0.2(a) and 
n = 0.4(b). 
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FIG. 4. The triplet pairing susceptibilities versus temper- 
ature for the triangular and t-t' lattices for n = 0.2(a) and 
n = 0.4(b). 




